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ABSTRACT 



We present new L-band spectra of 13 outbursting Be stars obtained with 
ISAAC at the ESO Paranal observatory. These stars can be classified in three 
groups depending on the presence/absence of emission lines and the strength of 
Bra and Pf7 emission lines relative to those of Humphreys lines from transitions 
6-14 to the end of the series. These groups are representative of circumstellar 
envelopes with different optical depths. For the group showing Bra and Pi'y 
lines stronger than Humphreys lines, the Humphreys decrement roughly follow 
the Menzel case-B for optically thin conditions. For the group showing compa- 
rable Bra, Pf7 and Humphreys emission line strengths, the Humphreys decre- 
ments moves from an optically thin to an optically thick regime at a transition 
wavelength which is characteristic for each star, but typically is located around 
3.65-3.75 iim (transitions 6-19 and 6-17). Higher order Humphreys lines probe 
optically thin inner regions even in the optically thicker envelopes. We find evi- 
dence of larger broadening in the infrared emission lines compared with optical 
lines, probably reflecting larger vertical velocity flelds near the star. The exis- 
tence of the aforementioned groups is in principle consistent with the proposed 
description by de Wit et al. (2006) for Be star outbursts in terms of the ejection of 
an optically thick disk that expands and becomes optically thin before dissipation 
into the interstellar medium. Time resolved L-band spectroscopy sampling the 
outburst cycle promises to be an unique tool for testing Be star disk evolution. 

Subject headings: stars: emission-line. Be - stars: mass-loss - stars: evolution - 
stars: activity 
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Introduction 



Be stars are rapidly rotati ng dwarf or giant B 



shown emissions in the Ha line ( jJaschek fc Jaschek 



type objects that show or have once 



that originate in the free- 


ree and free- 


30und 


(Gehrz. Hackwell & Jones 


1974; 


Waters 


1986) 



19871 ). Classical Be stars have moderate 
30und emission from ionized 



19861 ). Interferometric studies 



have shown that this gas is concentrated towards the equatorial plane forming a dense 
disk-like envelope extending extending up to ~ 10 stellar radii from the stellar surface 
(Grundstrom & Gies 2006, Quirrenbach et al. 1994, Stee et al. 1995). The IR region 
is domi nateci by broad and bright emission lines arising from high leyels of hydrogen 



atoms ( IBriot 



1981 



Andrillat. Jaschek. fc Jaschek 



1988 



Lenorzer et al. 



2002al ). He I, MgH 



and NaT emission lines in the K-band have also been reported (Clark and Steele 2000). 
he IR line optical depths and line flux ratios display a large variation fr om star to star 



Persson fc McGregor 



19851 ) and do not correlate with the spectral type (ILenorzer et al. 
2002al ). Be stars are also intrinsically variable, some of them show mild periodic or irregular 
photometric variability (Mennickent, Vogt & Sterken 1994, Sterken, Vogt & Mennickent 
1996), whereas others show sudden brightenings usually attributed to mass ejections from 
the surface of the stars, which can occur discretely over a range of timescales (Hubert, 
Floquet & Zorec 2000, Mennickent et al. 2002, de Wit et al. 2006). These outbursts 
probably induce variability in the opacity, the size and the geometry of the circumstellar 
envelope. Thus our aim is to explore the physical properties of envelopes of outbursting Be 
stars. To our purpose, we selected for an infrared spectroscopic study 6 Galactic Be stars 



showing long-lived outbursts (duration several hundred days) and 7 showing shor 



Hubert. Floquet fc Zored (120001 ). 



ived 



outbursts (duration days or tens of days) from the list of 
The stars were selected spanning a wide range of projected rotational velocities and most 
of them have been rarely studied spectroscopically. We hope contributing to the knowledge 
of the L-band spectral region in Be stars and its relation with the circumstellar envelope. 
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which have been hitherto scarcely studied. 

2. Observations and data reduction 

L-band spectra were obtained with the VLT Infrared Spectrometer and Array Camera 
(ISAAC, Moorwood et al. 1998) at the ESO Cerro Paranal Observatory in service mode 
during the nights of may 28-30, July 26, august 7-8 and September 3-4, 2003. The mode 
long slit low resolution spectroscopy was selected along with a central wavelength of 3.5 
fim. The pixel scale was 0.146 "/pixel. Two different setups were used, a narrow slit of 
0.3" and other of 2", providing resolving power of 1200 and 180, respectively. As our 
targets are bright, we had to use very short exposure times (0.3-2s) so that the detector 
does not become saturate; this had no effect on the total number of counts recorded, since 
the exposure time is controlled by the detector and no shutter geometrical bias is involved 
during the photon acquisition. Images were reduced with the ISAAC pipeline. The spectra 
were telluric corrected with the aid of early G-type telluric standards observed during the 
run at similar airmasses that science objects using the procedure described in Maiolino, 
Rieke & Rieke (1996). We built our telluric templates by dividing the telluric spectra 
by a synthetic solar-type spectrum interpolated at the same resolution and wavelength 
range. Then, we used the IRAlj^ telluric task to remove telluric absorption lines from the 
science objects. The telluric bands were successfully removed from our spectra, except in 
some cases in the wavelengths short-ward of 3.4 /im, characterized by heavy and variable 
atmospheric absorption. Nevertheless, it has minor importance since the spectral lines used 

^IRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Table 1: List of Be stars for which L-band spectra were obtained with ISAAC. A quahty 
flag (Q) is given for every spectrum: G for good and N for noisy 



Object N slit exptime airmass seeing MJD Q 

!"] W ["] 



V1448 Aql 


20 


0, 


,3 


1 


.85 


1 


.3 


0.6 


52789, 


.3411- 


52789, 


.3726 


G 


V1448 Aql 


6 


2, 


.0 


0, 


.28 


1 


.2 


0.7 


52789, 


.2966- 


52789, 


.3046 


N 


V618 Car 


4 


0, 


.3 


1 


.85 


1 


.5 


0.7 


52789. 


.0671- 


■52789, 


.0914 


G 


V618 Car 


2 


2, 


.0 


0, 


.28 


1 


.7 


0.8 


62789, 


.1062- 


■52789, 


.1078 


N 


V767 Cen 


2 


2, 


.0 


0, 


.28 


1 


.7 


0.7 


52860, 


.0731- 


■52860, 


.0746 


G 


V817 Con 


2 


0, 


,3 


1 


.85 


1 


,7 


0.8 


52847, 


,0313- 


52847 


.0329 


G 


V817 Cen 


2 


2, 


,0 


0, 


.28 


2, 


.0 


0.9 


52847, 


,0525- 


■52847, 


.0540 


G 


fx Cen 


2 


2, 


.0 


0, 


.28 


1 


.7 


1.0 


52859, 


.0684- 


■52859, 


.0700 


G 


fx Cen 


2 


0, 


,3 


0, 


.28 


1 


.6 


0.9 


62869, 


.0622- 


■52859, 


.0639 


G 


CP Cir 


4 


0, 


.3 


1 


.85 


1 


.6 


0.7 


62788, 


.2476- 


■52788, 


.2526 


G 


CP Cir 


2 


2 


,0 





.28 


1 


.4 


0.8 


52788, 


,2127- 


■52788, 


.2303 


N 


CP Cir 


2 


2, 


.0 


0, 


.28 


1 


.6 


0.6 


62788, 


.2643- 


■52788, 


.2659 


N 


KV Mus 


20 


0, 


.3 


1 


.85 


1 


.9 


0.7-1.0 


52868, 


.9961- 


■52859, 


.0246 


N 


OZ Nor 


2 


2 


,0 





.28 


1 


.1 


0.9 


52789 


,2050- 


■52789 


.2066 


G 


OZ Nor 


2 


0, 


,3 


1 


.85 


1 


.8 


0.9 


52790, 


,3767- 


52790, 


.3784 


G 


V467 Sot 


20 


0, 


.3 


1 


.85 


1 


.1 


0.7 


62788, 


.3288- 


■52788, 


.3603 


G 


V467 Set 


8 


2, 


.0 


1 


.85 


1 


.2 


1.1 


62789, 


.2186- 


■52789 


.2297 


N 


V341 Sge 


4 


0, 


.3 


1 


.85 


1 


.5 


0.4 


62789, 


.3886- 


■52789, 


.3918 


G 


V341 Sge 


2 


2 


,0 





.28 


1 


.4 


0.7 


52789, 


,3138- 


■52789 


.3155 


G 


V4024 Sgr 


2 


2, 


.0 


0, 


.28 


1 


.0 


0.7 


62789, 


,2791- 


■52789, 


.2806 


G 


V4024 Sgr 


2 


0, 


.3 


1 


.85 


1 


.8 


0.8 


62790, 


.3921- 


■52790, 


.3937 


G 


V1150 Tau 


16 


0, 


,3 


1 


.85 


1 


.2 


1.0-1.4 


52886 


4()()()- 


52886 


.4232 


G 


V1150 Tau 


6 


2, 


,0 


0, 


.28 


1 


.5 


0.8 


52887, 


,3249- 


■52887, 


.3324 


G 


V396 Vul 


2 


2, 


.0 


0, 


.28 


1 


.5 


0.7 


62789, 


.3268- 


■52789, 


.3274 


G 


V396 Vul 


2 


0, 


.3 


1 


.85 


1 


.7 


0.8 


62790, 


.4069- 


■52790, 


.4086 


G 
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in this work are mostly located out of this region. The spectra taken with the wide slit 
were flux calibrated with the aid of the standard star BS5471 (spectral type B3V), whose L 
magnitude is known. Those spectra taken with the narrow slit were continuum normalized. 
The observing log given in Table 1 indicates the number of spectra per object, single spectra 
exposure times and some additional observational parameters. 



3. Results 

The L-band IR spectral region, between 3.0 fim and 4.1 /im, displays Bra, Pf7, Pf6, 
Pie and most of Humphreys' lines from transition 6—14 to the end of the series. Weak 
Hel 4.038-4.041 /im emission is observed in /i Cen, OZ Nor, V4024 Sgr and V1150 Tau. 
Although the blue region was discarted by not providing a confident telluric substraction, 
we do observe Pf5 and Pfe emission in OZ Car, V817 Cen and V1150 Tau. The emission 
line spectra look similar to those of Be stars with irregular photometric variability, like 7 
Cas (Hony et al. 2000). According to the intensity of the emission lines, the sample of 13 
Be stars can be classified in three groups: those with Bra and Pf7 emissions equally intense 
as Humphreys' lines (Group I, gathered in Fig. 1); those with Bra and Pf7 emissions more 
intense than Humphreys' lines (Group II, see Fig. 2), and those where no line emission 
is detected (Group HI, Fig. 3). In this scheme, V1150 Tau could be a transition object 
between Group I and II. All stars show blue continua, and fi Cen is notable for showing a 
break in the continuum slope around A = 3.9 /im, which we attribute to an instrumental 
effect that should not affect our measurements of line intensity relative to the continuum 
and emission line widths. 

L magnitudes were computed by convolution of th e spectra obtaiii e d wit h the slit of 



2" with the transmission curve of the L filter given by 



Bessell fc Brett 



(jl988|). As part of 



the spectrum of fi Cen was corrupted and for KV Mus we had only one spectrum taken with 
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Fig. 1. — IR wide slit spectra of stars of Group I. 
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wavelenght (|j.m) wavelength (|a.m) 



Fig. 2.- 



IR wide slit spectra of stars of Group II and V 1150 Tau. 
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Fig. 3. — IR wide slit spectra of stars of Group III and the B3V star BS5471. 
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the narrow slit (0.3") which shows no emission hnes, no magnitude determinations were 
possible for these stars. In Table 2 we list the L magnitudes, the fundamental parameters of 
the program sta rs and their classification in the af orementioned groups. The V si n i value s 



were taken from 



f l200lh and 



Hubert 1 



Fremat et al 



(120061) 



loauet fc Zored fl2000h . Iciebocki fc Stawikowskil (|2000|) 



Yudin 



Some of our stars were too bright to det ermine the outbu rst stage (outburst/quiescence) 



from reliable ASAS-3 V-band light curves (jPojmahski 



200ll ) at the epoch of our L-band 



spectroscopy. Others were not included in such a catalog. The exception was OZ Nor, that 
was observed near maximum at the time of the L band spectroscopy. The Ha spectrum of 
V395 Vul taken in may 23, 2003 and of V4024 Sgr taken in may 29, 2003, suggest that our 
IR observations for these stars were obtained during minimum and maximum Ha emission, 
respectively ( http : / / astrosurf . com /bull /becat / ) . 



Line fluxes, full width at half maximum [FWHM) and equivalent widths, measured in 
the narrow slit spectra using the splot IRAF task, vary strongly from star to star. In Table 
3 we list these parameters for the stars showing emission lines (groups I and II). 
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Table 2: Spectrophotometric L magnitudes, fundamental parameters from given references 
and 2MASS color excess of the observed Be stars. The note indicates the classification group 
and the outburst character (1= long, s=short, see text). V sini is in kms~^. 



object L Teff(K)log^ S.T. V sini E{H - K) Note Rel 



V1448 Aql 7.66 


20000 


3.80 


B2 IVe 


243 


0.124 


III/s 


a 


V518 Car 


5.31 


18700 




B3-5Ve 


130 


0.043 


III/l 


c 


IJ, Cen 




23130 


4.04 


B2 IV- Ve 


155 




II/s 


d 


V767 Cen 


5.66 


23320 


3.95 


B31IIe 


70 


0.224 


I/s 


d 


V817 Cen 


5.25 


17490 


3.12 


B3IVe 


130 


0.178 


I/s 


b 


CP Cir 


8.00 


15 000 




B5IVe 




0.005 


III/l 


c 


KV Mus 




18 700 




B3Ve 




0.057 


III/l 


c 


OZ Nor 


5.64 


20300 




B2Ve 


86 


0.259 


I/l 


c 


V457 Set 


7.82 


23000 


3.80 


B1.5IVe 


202 


0.306 


III/s 


a 


V341 Sge 


7.58 


22 000 




B2.5Ve 


115 


0.105 


II/s 


c 


V4024 Sgr 


4.85 


18 940 


3.49 


B2Ve 


120 


0.113 


I/l 


d 


VI 150 Tau 7.25 


22 000 




B2.5Ve 


340 


0.323 


I/l 


c 


V395 Vul 


5.36 


17790 


3.93 


B2.5Ve 


230 


0.134 


II/s 


d 



(a) and (b) from Fremat et al. (2005,2006) 

(c) from Wright et al. (2003) ;(d) from Zorec et al. (2005) 
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Table 3: Spectroscopic measurement. We give the observed line fluxes / (± 10%, in ergs 
cm~^s~"^), equivalent widths EW (± 10%, in A) and full width at half maximum FWHM 
(± 40 kms~^). Ao defined in the text is in microns. 



object 


/-Bra 


fpfy 


fHuH 


-EWBr-c 


-EWpf^ 


-EWhuU 


FWHMBra 


FWHMpf^ 


FWHMH^,14, Ao 


/X Cen 


B.46E- 


■ 11 


2.32E- 


■ 11 


1.03E- 


■ 11 


82.8 


25.8 


16.0 


512 


370 


447 


> 4.0 


V767 Cen 


9.73E- 


■ 12 


1.17E- 


■ 11 


9.26B- 


■ 12 


28.3 


21.5 


24.7 








3.66 


V817 Cen 


1.62E- 


■ 11 


1.26E- 


■ 11 


1.23E- 


■ 11 


24.9 


17.2 


22.7 


277 


399 


367 


3.66 


OZ Nor 


1.2SE- 


■ 11 


9^50E- 


■ 12 


9.38E- 


■ 12 


35^0 


16.5 


23.3 


238 


288 


267 


3.75 


V341 Sgc 


3.41E- 


■ 12 


1.47E- 


■12 


6^8()E- 


13 


57.0 


17.7 


11.6 


136 


14,-) 


128 


> 4.0 


V4024 Sgr 


1.99E- 


■ 11 


2.23E- 


■11 


1.84E- 


■11 


23^3 


15^9 


19.8 


356 


327 


326 


3.75 


V1150 Tau 


B.91E- 


■ 12 


B^28E- 


■ 12 


4.04E- 


■ 12 


72.6 


44.2 


40.0 


350 


556 


567 


3.76 


V396 Vul 


1.06E- 


■ 11 


7^88E- 


■ 12 


6^68E- 


■ 12 


24.3 


14.3 


9.9 


638 


554 


484 


> 4.0 



13 



4. 



Discussion 



4.1. 



Kinematical issues 



Our line width accuracy is limited by the instrumental resolution, the signal to noise 
ratio of the spectra and the determination of the continuum level. We estimate the error 
of the full width at half maximum (FWHM) of the emission lines (corrected by the 
instrumental resolution and measured in the narrow slit spectra) to be ~ 40 km/s. On 
the other hand, the V sin i values given in the literature are usually accurate to less than 
~ 10%, although for fast rotators the values could be underestimated by larger amounts 
(Townsend, Owocki & Howarth 2004). Keeping in mind these caveats, we observe that 
FWHM roughly correlates with V sin i and also is larger for higher order transitions in the 
Humphreys series (Figs. 4 and 5). We found correlations of the form (Table 4): 



We considered only correlations with R > 0.60 and at most one rejected point. These 
correlations suggest that rotational broadening is the main line broadening mechanism for 
these lines and point to a rotationally supported, probably disk-like envelope as the source 
of infrared line emission. Similar kinematical insights have been derived from other infrared 
and optical spectroscopic studies of Be stars and they have been interpreted in terms of a 
disk-like geometry for Be star envelopes (e.g. Sellgren & Smith 1992, Hanuschik 1996, Clark 
& Steele 2000, Hony et al. 2000). In this view, higher order hues probe inner disk regions, 
with larger rotational velocities. The fact that the FWHM is larger than 2 V sin i could 
indicate additional sources of broadening like turbulence, macroscopic velocity fields or 
electron scattering (although we do not observe prominent electron scattering wings in the 



FWHM ^AxVsini + B. 



(1) 



-14- 



Table 4: Results of the FWHM — V sin i fits, it! is the correlation coefficient and N is the 
number of stars. 



Line 


A (km/s) 


B R a 


(km/s) 


N 


rejected 


source 


Pf7 


1.12±0.25 


218±49 0.83 


52 


6 


1 


a 


Hul4 


1.12±0.20 


209±39 0.89 


41 


6 


1 


a 


Hul5 


1.62±0.55 


163±103 0.63 


120 


7 





a 


Hul6 


2.04±0.33 


48±61 0.89 


71 


7 





a 


Hu21 


0.68±0.09 


363±18 0.95 


19 


5 





a 


Ha 


1.23 


70 




115 





b 


H/3+FeII 


1.2 


30 




90 





b 


Br7 


0.759 


149 




39 





c 


He I 2.058 


0.643 


247 




19 





c 



(a) this paper, (b) Hanuschik 1989, (c) Clarke & Steele 1999 
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Fig. 4. — FWHM versus V sini for the program stars for lines Pf7, Huu and Hu2i. For 
each panel, the upper dashed line is the best linear fit, the solid line represents the case 
FWHM — 2 V sini and the dotted horizontal hue indicates the spectral resolution (250 
km/s). 
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3,0 r- 

2,8 - 

2,6 - 

J, 2,4 - 

2,2 - 

I 1,8 - 

LL 

w 

^ 1,2 - 
X 

g 1,0 - 

0,8 - 




A 

_1_ 



— ■ — 


VII 50 Tau 


---0-- 


V341 Sge 


.... A ■ 


V395 Vul 


— V- 


|i Cen 




V817Cen 


M 


OZ Nor 





V4024 Sgr 



3,5 3,6 3,7 3,8 

wavelength (|xm) 




Fig. 5. — The normalized full width at half maximum versus wavelength for the Humphreys 
emission lines. 
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lines). We note that the B coefficient defined in equation (1) is the expected line broadening 
for a star seen pole-on {V sini = km/s). If planar Keplerian motions dominate the disk 
kinematics, then we should expect this number to be similar to the thermal broadening 
in the disk {Bth ~ 13 km/s for 10.000 K hydrogen gas). However, we observe all hues in 
Table 4 with B » Bth- Since IR hues, especially those of the Humphreys series, probe 
the inner disk region, it is possible that larger turbulent motions and eventually departures 
of the disk geometry in the inner disk explain why the B coefficients of infrared lines are 
larger than those of Balmer and Fell hues (Table 4). As we don't observe Stark broadened 
emission wings, we discard pressure effects as the cause for the large observed B. We failed 
to estimate the electron density of the envelope (and so an estimate of the pressure effect) 
using the Inglis- Teller formula since the quality of the data did not allow us to detect the 
wavelength of confluence of the Humphreys series with the continuum. The broadening 
effect mentioned, plus our rather small spectral resolution [R 1200), probably explains 
why the typical hallmarks for rotationally supported disks, viz. doubly peaked emission 
lines, are not observed in our spectra. 

In the disk model of Be stars the peak separation AA„ of the nth emission line measures 
the velocity Vn near the outer disk (Hirata & Kogure 1984): 

c^^VnSini (2) 

In the following we assume a disk rotational law given by: 



v(r) = v^(^) ^ r>r^ 



(3) 
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where is the equatorial stellar velocity and j is equal to 0.5 for a Keplerian disk and 1.0 
for continuum mass loss with conservation of the angular momentum. The extension of the 
disk which corresponds to the nth emission line is: 



and the extension relative to the Hul4 forming region is: 



Where we have assumed a linear relation between vhuu and FWHMhuu (Hanuschik 1996 
and references therein). The fact that the ratio [FWHMHun/ FWHMhuu) reaches values 
up to 2.7 around Hu24 (Fig. 5), implies that the relative disk extension rHui4:lfHu2i for OZ 
Nor and V341 Sge equals 3 (j=l) or 7 (j=0.5). 



4.2. Diagnostics for the envelope optical depth 



Lenorzer. de Koter. fc WatersI (l2002bl ) have suggested that the study of log (Hul4/BrQ;) 
versus log (Hul4/Pf7) provides a simple way to investigate the nature of the circumstellar 
material around Be stars and related objects. They found that line ratios close to the 
Menzel case B recombination theory correspond to optically thin lines formed mainly in 
isothermal stellar winds and those close to unity are optically thick lines arising from a 
disk-like structure. Particularly, in an optically thick medium, the line flux ratios become 
independent on the mass absorption coefficient and the line flux turns out to be dominated 
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by the size of the emitting surface. We show in Fig. 6 our decrements compared with those 
taken from the hterature or from theoretical models. It is clear that Group II stars are 
distributed in a region of moderate or small optical depth, whereas those of Group I are in 
the extreme of the distribution, corresponding to cases of optically thick envelopes. Group 
III stars, without emission lines, probably correspond to stars that have lost most of their 
envelopes. We note that not all of these stars show Bra and Pf7 absorption lines, suggesting 
the filling by emission in these lines and consecuently the presence of an incipient envelope. 
We find that our empirical groups trace the optical depth of the circumstellar envelope. 

We note the change in position of V395 Vul (= 12 Vul) in two epochs. As said in 
Section 3, this corresponds to a weakening of the Ha emission line strength, which is 



consistent with a more trasluce nt envelope in may 2003. We a 



spectra of the stars reported by 



Lenorzer. de Koter. fc Waters! (l2002bl ) and published by 



so iny estigated the ISO 



Vandenbussche et al. (2002). We note that above log (Hul4/Bra) ~ -0.2 all the stars (7 
objects, including V1150 Tau) can be classified as Group I, and below that hypothetical 
line the stars (9 objects) are Group II. 12 Vul changes between Group I and II. This shows 
that transition objects are difficult to find, and that the classification in groups I, II and III 
is i n principle observationally supported. This classification has a quantitative support in 



the iLenorzer. de Koter. fc WatersI (l2002bl ) diagram, but can be done more rapidly by simple 
visual inspection of the spectrum, for classification or selection purposes. We note that Be 
stars showing long-lived outbursts are found only in the upper right part of the diagram 
whereas those showing short-lived outbursts are more widely distributed. Our qualitative 
conclusion about the optical depth condition in the envelopes of Group I & II stars is also 
supported by the emission line ratio Pf7/Hul6, which is a good discriminator between 
optically thin and thick conditions (Hummer & Storey 1987, Hamann & Simon 1987). We 
found this ratio ~ 1 for all Group I stars and much larger than unity for Group II stars, 
consistent with theoretical predictions for optically thick and thin envelopes (respectively) 
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with r= 10.000 K and 10^° cm-^ 

A quantitative interpretation of the Lenorzer et al. diagram was done by Jones et al. 
(2009), in terms of a disk with varying density and illuminated by a central star of given 
temperature Tg//. In this model the disk density profile is given by: 

p(i?,Z)=po(f)"e-(W (6) 



where R is the cylindrical, radial distance from the star rotation axis, Z is the perpendicular 
height above the equatorial plane and H the corresponding scale height, po is the density 
at the inner edge of the disk in the equatorial plane. Jones et al. (2009) explain the region 
occupied by Be stars in the Lenorzer et al. diagram by different values of po; ^ and T^ff. 
They note that the models correctly explain the observed decrements if n = 3.0 to 4.5 and 
Po— 7.5x 10~^^ to 1.0 X 10~^^. Their models do not reproduce, however, the extreme 
optically thick cases located at the upper right corner of the diagram (our Group I stars). 
The fact that in Group I stars Bra is comparable in strength to Humphreys lines suggests 
that the size of the emitting region is limited (by a putative companion for instance) and 
the gas density is high. However, the fact that we observe Bra line widths lower than 
in Humphreys lines for all Group I cases, argues against this view. Also, the Jones et 
al. (2009) disk models cannot reproduce these cases by simply truncating the disk since, 
as reported by these authors, all line fluxes fall off approximately at the same rate with 
increasing distance from the central star. Jones et al. (2009) claim that these stars must 
have disk geometries different than a simple power-law. In this sense, different geometries 
like power-laws with non constant n (Zorec et al 2007) could be explored. Our discovery 
of a significant number of Be stars in this region calls for additional theoretical work for 
explaining these cases. 
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4.3. Study of the EWs and Humphreys emission Une decrements 

We find a strong correlation between EW/ X and wavelength for a given series. This 
parameter increases with A and sometimes saturates at EW/X fa 6 x 10"^ (Fig. 7). This 
kind of behavior was observed in the prototype Be star 7 Cas in several infrared H I series 
and interpreted in terms of a decrease in the line source function in the outer parts of the 
circumstellar region (Hony et al. 2000). Here we demonstrate that this is an usual behavior 
of Group-I and Group-II stars, being the pattern only disrupted by the larger emission 
found in the Bra and Pf7 lines of Group-II stars. 

We have studied f{Hun)/f{Hul9) (n is the quantum number of the upper level) 
versus A for every star showing emission lines. We find that the Humphreys decrements 
follow well defined patterns (Fig. 8). For the stars of Group II, they increase with A, 
indicating optically thin conditions, as in the case of P Monocerotis A (Sellgren & Smith 
1992). On the contrary, for stars of Group I, the decrements change their behavior at 
a transition wavelength, Aq, from optically thin conditions (at shorter wavelengths) to 
optically thick conditions (at longer wavelengths). This transition is rather fast, and Aq, 
given in Table 3, apparently is a characteristic of every star. As in the case of EW/ A, the 
transition at Aq could reflect a change in the optical properties of the envelopes at a certain 
distance of the star. In this scheme, V1150 Tau could be classified as a Group-I star. We 
note that even for Be stars with optically thick envelopes, high order Humphreys lines 
probe optically thin inner regions. 

From Fig. 5 and 7 we deduce a rough anticorrelation between FWHM and EW/ X 
for the Humphreys lines of a given star, as usually happens for the Ha line in Be stars 
(Dachs et al. 1986). Consequently, this can be explained assuming that the EW scales with 
the disk size and that higher order lines are formed in smaller and inner disk-like regions 
rotating faster than the outer and bigger regions forming low order Humphreys lines. 
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Fig. 6. — Line ratio diagram for our objects (stars) and Be stars from Lenorzer et al. (2002a) 
(dots). The solid square indicates the loci for an optically thick black body and the big disk 
the range of ratios for Menzel case B recombination, including collisional de-excitation, for 
temperatures higher than 10^ K. Dashed lines illustrate the path of the models by Jones et 
al. (2009) for a 12 Mq B2 star with luminosity 7000 L©, T^ff= 20.000 K, radius 7.0 i?©, log 
g= 4.0, inclination angle 35 degree and surrounded by a disk of solar composition, exponent 
n= 3.5 and increasing density po from point 1 to 2. The model in the upper position has 
Po= 10^^^ g cm~^. Above log (Hul4/Bra) -0.2 all the stars can be classified as Group I 
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4.4. 2MASS photometry and continuum slope 

For each star we have calculated E(H — K) — [H — K)obs — {H — K)q, where the first 
color is from 2MASS (Skrutskie et al. 2006) and the last one is given by Koornneef (1983) 
for a star of similar spectral type and luminosity class (Table 2). We have not corrected 
them by interstellar extinction, which is expected to be small in infrared wavelengths, 
especially for our bright and relatively close stars. Keeping in mind that variability could 
affect the comparison of non-simultaneous data, we find no correlation between color excess, 
outburst character and group membership. 

We measured the flux ratio in regions almost depleted of emission lines: 



flux ratio = (7) 
J3.93 ^>^'^^ 



where S\ is the spectral flux density. This ratio is bluer for hotter stars and for a given 
temperature is bluer for stars without emission lines (i.e. Group 3 stars. Fig. 9). No 
correlation is observed between continuum slope and emission line strength. These findings 
possibly point to the importance of the stellar flux in the infrared continuum emission; 
hotter stars produce bluer continuum and when the disk becomes developed (Group I-II 
stars) this continuum becomes redder, probably due to the contribution of free-free emission 
and partial stellar obscuration from the optically thicker disk. The position of V457 Set 
is notable. This star shows the bluest continuum, and no developed emission. However, 
the absence of Bra and Pf7 absorptions indicate that residual emission is filling these 
lines. Incidentally, V 1448 Aql also show these characteristics; a very blue continuum and 
filled absorption lines. We speculate that these stars could be in a special position of their 
eruptive cycles, such as in the process of ejecting a hot optically thick disk. The disk could 
be later dissipated as an optically thin ring into the interstellar medium, as proposed in the 
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model of de Wit et al. (2006), causing the stars' position to move down in the diagram. The 
bluer color during the outburst ascending branch are predicted by this model, and observed 
in Be outbursting stars (de Wit et al. 2006). The stars should move up and down in the 
diagram of Fig. 9 during the outburst rising and decay. We note that Be star outbursts 
seem to be of larger amplitude in red bands (Mennickent et al. 2002) so this effect in our 
studied region should be significant. Further time resolved spectroscopy sampling the whole 
eruptive cycle is needed to test this conjecture for Be stars. 



4.5. Dust eiround outbursting Be stars? 

We note the absence of dust spectral features at the L band, especially the Polycychc 
Aromatic Hydrocarbons (PAH) emission feature at 3.3 //m and the nano-diamond features 
at 3.43 and 3.52 /xm, that have been detected in some Herbig Ae/Be stars. These features 
have been observed in pre-main sequence stars of spectral type B9 and later but for earlier 
spectral types the 3.3 //m band is weak or absent and the other bands are even weaker 
(Habart et al. 2004, Acke & van den Ancker 2006). As our stars are mostly early B-type, 
we cannot establish the absence of dusty envelopes from the absence of these key spectral 
features. However, the infrared color excesses E{H — K) of our targets hsted in Table 2 
and the 2MASS colors, J -K between -0.08 and +0.34 and - X between -0.29 and 0.30, 
arc typical for Be stars and not so large as in most hot stars surrounded by dust (e.g. Fig. 
11 in Mathew, Subramaniam & Bhatt 2008). This suggests the absence of significant dust 
in the envelopes of these ourbursting Be stars. The IRAS colors of our targets provide the 
same insight. This finding is consistent with the classification of our targets as "canonical" 
Be stars. Our targets probably are not surrounded by massive cool envelopes as in the case 
of Herbig Ae/Be stars. 
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5. Conclusions 

We have provided a view of the L-band spectra of a selected sample of outbursting Be 
stars. These spectra show no evidence of dust and are not different from those reported 
spectra of Be stars showing only irregular photometric variability. The observed L-band 
spectra of 13 outbursting Be stars can be categorized in three broad groups reflecting 
the optical depth conditions in the Be star envelope. Based on the relative intensity of 
Humphreys, Bra and Pf7 emission lines, a rapid visual inspection of the spectra indicates 
the envelopes optical depth. In addition, the Humphreys decrements, and the parameter 
Ao defined in Section 4.3, can be used as a diagnostic tool for the optical depth of the 
circumstellar envelopes. We find that higher order Humphreys lines probe optically 
thin inner regions even in the case of optically thick envelopes. In addition, the large 
broadening observed in the IR lines probably reflect vertical velocity fields near the star. 
Some power-law disk models describing the infrared emission line properties fail to explain 
the cases of optically thick envelopes (Group 1 stars). We expect that our discovery of a 
large number of these stars motivates further theoretical work in this area. Our data does 
not allow us to test a possible correlation between the outburst stage and the spectral 
appearance, but the fact that the stars were observed at random outburst phases, the 
changes observed in two epochs in V395 Vul (= 12 Vul) and the blue continuums of 
V1448 Aql and V457 Set, suggest that all outbursting stars observed in this project could 
pass during their cycles through Groups I, II and HI. The existence of these groups is in 
principle consistent with the proposed outburst description by dc Wit et al. (2006) in 
terms of the ejection of an optically thick disk that expands and becomes optically thin 
before dissipation into the interstellar medium. Accordingly, Aq, and the whole spectral 
appareance should change in a significant way during the entire outburst cycle, following the 
development of the circumstellar envelope. Variability of Be stars in the diagram of Fig. 6 
along the diagonal was already suggested by Lenorzer et al. (2002b) due to the transient 
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nature of Be star disks. We further suggest that the Be star outburst mechanism and disk 
evolution could be tested with time resolved infrared spectroscopy, using for instance the 
spectral diagnostics worked in this paper. Further time resolved L-band spectroscopy of 
selected targets along with photometry spanning several eruptive cycles are needed to test 
our prediction for outbursting Be stars. 
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